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I. 


INTRODUCTION 


Modeling  the  Ignition  and  early  combustion  of  a solid  propellant 
granular  bed  has  typically  assumed  that,  once  Ignited,  the  regression 
rate  depends  only  on  the  Instantaneous  gas  pressure  (usually  r*aP”). 

For  lumped  parameter  gun  codes,  like  Baer-Fr ankle It  has  been  the 
practice  to  adjust  the  two  parameters  (a,  n)  In  the  burning  rate  ex- 
pression to  Induce  an  agreement  between  measured  gun  pressures  and 
computed  pressures.  Such  an  adjustment  process  permits  the  burning 
rate  to  compensate  to  some  degree  for  other  errors  In  input  data.  Al- 
though the  adjustments  do  permit  parametric  studies  of  gun  performance, 
the  disagreement  with  Independent  burning  rate  measurements  prohibit 
the  usefulness  of  the  burning  rate  formula  for  any  other  gun. 

Within  the  past  five  years,  there  have  been  several  one-dlmenslonal 
models  of  the  flame  spreading  and  combustion  of  the  granular  propellant 
bed.  All  these  models  use  the  same  power  law  burning  rate  formulation 
as  the  lumped  parameter  models  (r^aP'^) . Until  recently  no  serious  con- 
sideration has  been  given  to  using  other  laws  for  burning  rates  pro- 
bably because  either  no  definitive  solutions  are  available  for  the 
burning  rate  and/or  the  approximate  models  require  too  much  additional 
computing.  To  date  there  have  not  even  been  approximations  used  to 
find  a model's  sensitivity  to  potential  transient  burning  effects. 

Considerable  work  has  been  directed  to  transient  combustion  of 
solid  propellants  in  rocket  motors.  Emphasis  has  been  on  linear  or 
nearly  linear  perturbations  in  chamber  pressure  oscillation  or  in  ex- 
tinguishment by  depressurization.  Kooker  and  Zinn^  did  a numerical  In- 
vestigation of  the  transient  burning  coupled  to  a solid  rocket  motor 
chamber  by  a nonlinear  method  and  predicted  some  sharp  burning  rate 
transients  with  only  small  amplitude  pressure  oscillations.  Recent 
work  by  Caveny  et  al^  attempted  to  model  transient  burning  In  gun  con- 
ditions by  using  an  approach  due  to  Zeldovich  for  transient  heat  feed- 
back to  the  propellant  surface.  Nonlinear  models  of  the  combustion 


P.  G.  Baer  and  J.  M.  Frankie^  "The  Simulation  of  Interior  Ballistia 
Performance  of  Guns  by  Digital  Computer  Program",  BEL  Report  No. 

1183,  Deaemher  1962  (AD  if 299980). 

D.  E.  Kooker  and  B.  T.  Zinn,  "Numerical  Investigation  of  Non~Linear 
Axial  Instabilities  in  Solid  Rooket  Motors",  BRL  Contract  Report  No. 
141,  March  1974  (AD  §776954). 

L.  H.  Caveny,  M.  Summerfield,  and  C.  W.  Nelson,  "Ignition  Transients 
and  Pressurization  in  Closed  Chambers",  11th  JANNAF  Combustion  Meeting, 
Pasadena,  CA,  CPIA  Pub.  261,  September  1974. 
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were  developed  by  Krler  et  al  and  by  Levine  and  Cullck  . When  Kooker 
and  Nelson^  compared  the  predictions  of  the  KTSS,  Levine,  and  Kooker 
models  In  a typical  gun  pressure  environment,  they  predicted  that  three 
models  had  the  same  qualitative  response  of  a single  sharp  transient  at 
a low  pressure  and  an  asymptotic  approach  to  the  power  law  represen- 
tation as  pressure  Increased. 


It  Is  the  purpose  of  this  brief  report  to  compare  the  burning 
rate  predictions  of  three  different  types  of  models  for  an  Imposed 
pressure  history  found  in  a typical  Army  gun. 


II.  MODELS 


dt 


Dependence  Models 


Early  approximations  for  the  dynamic  response  were  proposed  by 
Paul  et  al^,  Parker  and  Summerfleld°  and  vonElbe^.  With  some  simplified 
analyses  these  researchers  all  arrived  at  a form  for  instantaneous 
burning  rate: 


.no 


2 dt 


]. 


H.  Krier,  J,  S.  Tien,  W,  A,  Sirignano,  and  M,  Surmevfield,  "Nonsteady 
Burning  Phenomena  of  Solid  Propellants:  Theory  and  Experiments"  AIAAJ 

Vol.  6,  No.  2,  pp  278-285,  February  1968. 

J,  N.  Levine  and  F.  E.  C.  Culiak,  "Nonlinear  Analysis  of  Solid  Rocket 
Combustion  Instability",  AFRPL  TR-74-45,  October  1974. 

D.  E.  Kooker  and  C.  (/.  Nelson,  "Numerical  Solution  of  Three  Solid 
Propellant  Combustion  Models  During  a Gun  Pressure  Transient",  12th 
JANNAF  Combustion  Meeting,  Newport,  RI,  August  1975. 

B.  E.  Paul,  R.  L.  Levine,  and  L.  Y.  Fong,  "A  Ballistic  Explanation  of 
the  Ignition  Pressure  Peak",  AIAA  Solid  Propellant  Rocket  Conference, 
Palo  Alto,  CA,  AIAA  Paper  64-121, 

K.  Parker  and  M.  Surmerfield,  "The  Ignition  Transient  in  Solid  Pro- 
pellant Rocket  Motors",  Princeton  University,  Dpt.  AMS  Report  769, 
1966. 

G.  vonElbe,  "Solid  Propellant  Ignition  and  Response  of  Combustion  to 
Pressure  Transients,  AIAA  Paper  66-668,  1966. 
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where 


r 


a 

n 

P 


dt 


steady  state  burning  rate 
a constant  coefficient 
thermal  dlffuslvlty 
burning  rate  exponent 
Instantaneous  pressure 

time  derivative  of  pressure. 


Krler  briefly  reviewed  the  derivations  and  their  limitations  to  small 
changes  from  steady  state  operation.  The  coefficient  tj;  is  different 
for  each  derivation  (1/2  for  Parker  and  Summerfleld;  1 for  Paul  et  al; 

2 for  vonElbe) . Krler  then  relaxed  the  assumption  of  constant  surface 
temperature  and  used  the  KISS  model  of  dynamic  heat  feedback  to  derive 
a variable  expression  for  i|;: 


[1  _ -2_] 

^ _ m , 

[(p/Po)”^”(2+l/n)-2H] 

where  P reference  pressure 
m a constant 

H a non-dimensional  surface  heat  release. 

Quasi  Steady  Heat  Feedback  One  Dimensional  (QSHOD) 

Three  models  solve  the  partial  differential  energy  equation  in  the 
solid  propellant  using  a nonlinear  boundary  condition  at  the  propellant 
surface.  Kooker  and  Nelson  rederlved  the  governing  equations  in  a com- 
mon coordinate  system  and  compared  the  solutions  for  the  pressure  field 
in  a 105mm  tank  gun  as  is  the  purpose  of  this  report.  The  models  all 
use  a surface  boundary  condition  of  the  form 


dx 


+ G2/R  , 


where  ^ non  dimensional  thermal  gradient 
R non  dimensional  regression  rate, 

but  differ  in  the  form  of  the  functions  G^^  and  G2  and  R.  Table  I gives 
the  differences. 


ff.  Kviev,  "Solid  Propellant  Burning  Bate  During  a Pressure  Transient", 
Comb  Sai  Teah  Vol  5,  pp  69-73,  1972. 
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Table  I.  Heat  Feedback  Terms 


Model  ^1  S R 


KTSS 

H 

2n.  n/m 
P (p  -H) 

0 ” 
s 

Kooker/Zlnn 

C 

H+(0  -1)(1-  -£) 
s C 

s 

S SO  0 

A(0^-1) 

exp[  ] 

Levlne/Cullck 

c 

H+(0g-l)(l- 

s 

p^"{  f (-^ 1)  1 

^ 1-^lnP 

+(1-H)1 

A(0  -1) 
exp  [ ] 

where  0 non  dimensional  temperature 
c heat  capacity 
T^^  reference  surface  temperature 

A,E  non  dimensional  kinetic  constants. 

The  G]^R  term  represents  heat  from  the  exothermic  surface  reaction  and 
G2R  represents  heat  from  the  flame.  Kooker  and  Nelson  found  that  the 
Kooker/Zlnn  and  Levlne/Cullck  models  gave  practically  the  same  response 
and  that  the  KTSS  model  predicted  higher  overshoots  because  of  the  pre- 

C 

sence  of  the  (0  -1) (1-  term  In  the  surface  heat  release  term. 

s 

Zeldovlch  Model 


Caveny  et  al  used  the  Zeldovlch  approach  to  derive  a nonsteady 
heat  feedback  of  the  form 


dx 

o 


-[T  -ln(-r^)/a  +T  ] 
a s R p o 
o 


where  o 

P 


temperature  sensitivity. 


'i 


i 


f 
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The  advantage  gained  la  that  no  model  of  the  flame  need  be  specified. 
The  non-steady  heat  feedback  from  the  flame  at  R(Tq,P)  is  assumed  to 
be  the  same  as  the  feedback  at  RoC^oeq*^^  when  Tq-o  Is  the  Initial 
temperature  which  would  yield  steady  burning  at  R(P) • With  that 
boundary  condition, the  thermal  energy  equation  of  the  solid  propellant 
is  solved  numerically. 

The  Test  Case 


A reference  set  of  thermal  and  thermodynamic  properties  were 
assigned  to  a double  base  propellant.  Best  estimates  from  whatever 
experimental  data  were  available  were  used.  The  only  parameter  which 
is  observed  to  have  a dramatic  effect  on  the  results  is  the  value  of 
the  surface  heat  release  in  the  QSHOD  model  group.  Experimental  data 
from  Kubota  et  al^^  establish  that  the  value  can  range  easily  from  70 
to  100  cal/gm  which  corresponds  to  values  of  the  parameter  H from  0.58 
to  0.83. 

Temperature  sensitivity  for  the  propellant  was  measured  as  .DOAb/^K. 
To  be  consistent  with  the  other  parameters  including  H,  a rough  estimate 
of  the  temperature  sensitivity  may  be  obtained  from  Summerfield  et  al^2 
as  0.0054/°K.  One  source  of  the  difference  is  the  derivation's  assump- 
tion that  the  KTSS  flame  contribution  to  the  heat  feedback  is  independent 
of  initial  temperature.  The  difference  carries  no  immediate  impact  be- 
cause the  results  are  essentially  alike  for  either  value. 

The  imposed  pressure  profile  was  that  of  the  105mm  tank  gun 
starting  at  6.9  MPa  where  the  propellant  is  assumed  to  be  burning  at 
steady  state. 


III.  RESULTS  AND  DISCUSSION 

The  shapes  of  the  burning  rate  response  of  the  three  types  of 
models  are  distinctly  different.  Figure  1 plots  the  relative  burning 
rate  against  time  in  a pressure  field  which  rises  almost  linearly  at 
an  average  rate  of  about  10^  MPa/sec.  For  the  models  using  surface 
heat  release,  it  was  assumed  to  be  96  cal/gm. 


N.  Kubota,  T.  J.  Ohlemiller,  L,  H.  Caveny,  and  M.  Swmerfield,  "An 
Experirrtental  Study  of  the  Site  and  Mode  of  Action  of  Tlatonizevs 
in  Double  Base  Propellants" , AIAA  Paper  74-124,  January  1974. 

M.  Swmerfield,  L.  H.  Caveny,  R.  A.  Battista,  N,  Kubota,  Yu  A. 
Gostinsev , and  H.  Isoda,  "Theory  of  Dynamic  Extinguishment  of  Solid 
Propellants  with  Special  Reference  to  Nonsteady  Heat  Feedback  Law", 
Journal  of  Spacecraft  Rockets,  8_  13),  March  1971,  pp  251-258. 
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Two  values  of  ip  in  the  dp/dt  model  showed  different  results,  al- 
though the  same  general  shape.  The  Krler  version  reached  a peak  of 
12.4  when  the  vonElbe  version  reached  17.4.  The  vonElbe  curve  was  con- 
sistently higher  than  the  Krler  curve  throughout  the  time  examined. 

For  the  Krler  version,  the  coefficient  was  initially  about  1.5  and 
declined  steadily  to  0.7  at  0.5ms. 

The  QSHOD  models  all  had  the  same  shape  which  produced  a sharp 
peak  at  0.25ms.  Peaks  for  the  Levine/Culick  and  Kooker/Zinn  were  about 
the  same  at  6.1  but  the  KTSS  peak  was  much  higher  at  about  17  (not  shown 
in  Figure  1).  Kooker  and  Nelson  provide  a more  complete  comparison  for 
this  type  model. 

For  the  Zeldovlch  model  a much  lower  response  was  observed.  The 
peak  overshoot  of  only  2.1  occurred  at  0.4ms. 

In  the  computations,  burning  rate  excursions  did  not  affect  the 
imposed  pressure.  In  an  actual  chamber  the  two  would  be  coupled  be- 
cause the  faster  gas  generation  rate  would  increase  the  pressurization 
rate.  The  mass  of  gas  in  the  chamber  is  proportional  to  the  Integral 
/^rdt.^  To  the  extent  that  this  integral  exceeds  the  quasi-steady  inte- 
gral / r^dt,  the  chamber  pressurization  will  be  altered  by  the  burning 
rate  excursions.  A rough  estimate  of  the  effect  is  given  in  Table  II 
where  the  ratio  /rdt//r  dt  is  compared  at  different  times.  Translating 
the  ratios  into  real  chamber  effects  would  require  a coupling  to  the 
chamber  dynamics. 

Experimental  verification  of  the  magnitude  of  such  burning  rate 
excursion  has  not  been  made.  Closed  vessel  tests  of  granular  beds 
typically  discard  the  low  pressure  Information  because  of  uncertainties 
in  surface  area  and  ignition  spreading.  Sharp  spikes  have  been 
only  occasionally  reported  and  usually  Ignored  as  a test  anomally. 

An  inference  to  be  drawn  is  either  the  spikes  do  not  occur  or  the  time 
for  flame  spreading  is  of  the  same  order  as  the  time  for  the  burning 
rate  excursion.  Resolution  awaits  better  experimental  information. 


Table  li.  Gas 

Production 

Ratios 

0. 1ms 

0.25ms 

0.5ms 

0.8ms 

dp/dt  \p  = 2 

13.1 

7.10 

3.91 

2.47 

dp/dt  tp  = f (p) 

9.0 

4.75 

2.64 

1.79 

Levine/Culick 

0.97 

2.16 

1.81 

1.52 

Zeldovlch 

0.97 

1.21 

1.43 

1.40 

A serious  criticism  may  be  leveled  at  applying  the  QSHOD  models 
to  nitrocellulose  based  gun  propellants.  The  derivation  of  the  flame 
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contribution  to  the  heat  feedback  assumes  a uniform  distribution  of 
the  energy  release  in  the  flame  on  the  basis  that  composite  propellant 
flatues  are  diffusion  controlled.  Because  nitrocellulose  based  propellants 
are  avidly  believed  to  have  premixed,  klnetlcally  controlled,  thin 
flames,  the  models  must  be  reexamined  to  Judge  the  effect  of  changing 
the  flame  structure.  A quick  calculation  for  a flame  sheet  approximation 
showed  the  overshoots  to  be  more  like  the  Zeldovich  results.  The  true 
energy  release  distribution  lies  probably  somewhere  in  between  the  two 
extremes  with  a reasonable  conclusion  that  no  clear  case  can  be  made 
for  choosing  between  the  Zeldovich  and  QSHOD  models  for  a nitrocellulose 
based  propellant. 

The  implications  of  this  study  are  not  limited  to  the  particular 
problem  chosen.  The  gun  pressure  profile  selected,  monotonic  and 
nearly  linear,  is  not  strictly  limited  to  guns.  Two  real  cases  can  be 
conceived  where  the  results  are  applicable:  A gun  with  a composite 

propellant,  like  a nitramine  base,  or  any  other  vessel  like  a rocket, 
with  a composite  propellant  in  a rapid  pressurization.  The  message  is 
clear:  do  not  choose  a model  on  the  basis  that  they  all  give  the  same 

answer. 
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